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1 Introduction 


One of the main concerns in constructing a physically acceptable held theory is to avoid 
unphysical propagating modes, like ghosts. General Relativity (GR) is one example of 
a healthy classical theory of gravity which is based on a unique spacetime metric along 
with its Levi-Civita connection, and it is constructed out of second derivatives in the 
metric tensor. The simple form of Hilbert-Einstein action leads to held equations with no 
more than second derivatives, resulting in a gravitational theory without ghosts [1], which 
moreover, passes successfully with high accuracy all local observational tests both for weak 
and strong gravity [2]. 

Torsion is an important tensor besides curvature, which can be used for the construc¬ 
tion of gravitational theories. In this case, one can use as dynamical variables the vielbein 
and the connection, and the held equations arise by variation with respect to these two 
variables. Soon it was realized that imposing the teleparallelism condition, it is possible 
to express equivalently the standard Einstein gravity [3-10] (or even Gauss-Bonnet grav¬ 
ity [11]) by a different geometry described in terms of torsion. This can be achieved by 
employing instead of the Levi-Civita connection the Weitzenbdck one, which has vanishing 
curvature but non-vanishing antisymmetric piece and therefore non-vanishing torsion, and 
thus the only dynamical variable that remains is the vielbein. 

In a more general formulation of teleparallelism, an arbitrary connection of vanishing 
curvature is assumed, and everything is again expressed in terms of the torsion, but now 
both the vielbein and the connection are dynamical. In this formulation, the theory is both 
Lorentz and diffeomorphism invariant, while the former formulation lacks local Lorentz in¬ 
variance. In its teleparallel formulation, Einstein gravity action consists of the torsion 
scalar T, which is a specihc combination of quadratic torsion scalars and contains up to 
first-order vielbein derivatives. Then, the variation of this action gives exactly the Ein¬ 
stein gravitational field equations, and for this reason the theory is termed as Teleparallel 
Equivalent of General Relativity (TEGR) [3-10]. 

One interesting question is what kind of theories can be produced in which the tor¬ 
sion can give sizeable physical effects that can differentiate these theories from GR. One 
approach that was put forward, inspired from the f{R) modifications of GR, is to start 
from TEGR instead of GR, and construct f{T) extensions. This new class of gravitational 
modification do not coincide with f{R) and proves to have interesting local black hole 
solutions [12-21], and it gives variable cosmological models [22-33]. 

Another approach to probe the effects of torsion to GR, is to couple the torsion to 
scalar fields. Self-gravitating scalar helds were used to study the dynamics of local black 
hole solutions of Einstein equations. These studies had resulted to the well known no-hair 
theorems [34, 35] and to the way to evade them [36-40]. In cosmology the scalar fields 
have been also used to investigate the features of the early universe [41-43]. Apart from 
the simple minimal coupling case, the scalar field may be coupled non-minimally with 
curvature, since the appearance of terms of the form f((j))R is motivated by many reasons, 
such as the variability of the fundamental constants, the Kaluza-Klein compactification 
scheme, or the low-energy limit of superstring theory. 
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To capture the effects of curvature in the early stages of cosmological evolution or near 
the horizon of a black hole, couplings of scalar fields directly to curvature through terms of 
the form dyCp, were studied. However, these derivative coupling terms lead to held 

equations containing higher than second derivatives [44], and therefore ghosts. To remedy 
this pathology one has to couple the scalar helds to the Einstein tensor [45] and this term 

d^pducj) is part of the general Horndeski Lagrangian which is known to result [46] in 
a healthy theory with only second order held equations. These theories have interesting 
cosmological implications [45, 47-49] and they give also local black hole solutions [50-54]. 

In this work we will discuss hrst couplings of torsion to scalar helds of the form 
We will show that such coupling in four or higher dimensions, for a diagonal vielbein, does 
not give non-trivial spherically symmetric solutions^. We will then introduce a derivative 
coupling of a scalar held to torsion of the form Tg^^^d^pducj). This coupling, contrary to its 
curvature counterpart Rg^’-'d^j^pd^p-, has the advantage that it leads to second-order equa¬ 
tions of motion. We will show that such a coupling gives spherically symmetric solutions 
which have distinct features from their GR counterparts. Note that the above derivative 
coupling between torsion and the scalar held is not the only one that can be constructed, 
since other combinations of with quadratic (of even higher) powers of torsion could 

appear. 

Our work is organized as follows. In Section 2, after a brief review of scalar-tensor 
theories, we set up our theory of a scalar held coupled to torsion with a derivative coupling 
and present the held equations of the gravity-scalar system. In Section 3 we consider a 
spherically symmetric ansatz for the metric and we reduce the coupled held equations to a 
decoupled master equation from which the solutions are obtained. In Section 4 we discuss 
an exact wormhole-like solution, while in Section 5 we discuss large- distance solutions. 
Finally, Section 6 is devoted to our conclusions. 

2 Derivative coupling of a scalar field to torsion 

In this section, after a brief review of scalar-tensor theories, and in particular of a scalar 
held minimally coupled to gravity, we will introduce torsion, and we will set up the gravity- 
scalar held equations describing the torsion held coupled with a derivative coupling to a 
scalar held. 

In D dimensions the Einstein-Hilbert action with a scalar held minimally coupled to 
gravity is given by 

Smin = J d^xy/\i\R- J d^xy/\i\ dyp + v'^ . ( 2 . 1 ) 

The Ricci scalar R is constructed from the Christoffel connection Indices a,b,... will 

refer to tangent space, while /r, ... are coordinate ones. A cosmological constant A = 

±3Z“^, where I is the length of the (A)dS space, has been incorporated in the potential, so 

^The cosmological implications of these type of couplings were studied in [55-59] 

^On the contrary, in three dimensions spherically symmetric solutions were found [60]. 
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the total potential is V(</>) = -^+Vi{4>)- In four dimensions and for a negative cosmological 

Kd 

constant this theory has a local black hole solution known as MTZ solution [36] given by 


ds^ = B{r) 


F{r)dt^ + 


—^dr‘^ + 
F{r) 


( 2 . 2 ) 


where 


B{r) 


r(r + 2GAr/r) 
(r+GATAt)^ 




2 


(2.3) 


The constant /i can be identified with the mass of the black hole and «; 4 ^ = SvrG/v, with 
Gv the Newton’s constant. The line element da"^ refers to a two-dimensional manifold of 
negative constant curvature. The scalar held is given by 


0 = a/T-t;— arctanh ("——') , (2.4) 

V 47rGAr Vr + GAr/i/ 

and the potential I^i is found to be 

This is the simplest known hairy black hole solution of a scalar held minimally coupled 
to the curvature, which goes to zero at inhnity and it is regular on the horizon. Here we 
should stress an important property of this solution. If we switch off the scalar held we 
cannot get a black hole solution with hxed mass. In fact, there is only one integration 
constant /i, and for </>—)• 0 the geometry approaches a massless black hole. This means 
that for a given mass there are two branches of different black hole solutions, the one with 
non-trivial scalar held given by (2.2), and the vacuum solution (with = 0). 

In writing (2.1) we used the usual Christoffel connection which by construction has 
zero torsion, and thus, all the gravitational information is incorporated by the curvature 
(Riemann) tensor. In the teleparallel formulation an alternative connection is adopted, 
and the basic object is the torsion 


be ^ cb ^ be ^b ^ 


( 2 . 6 ) 


expressed in terms of the arbitrary vielbein Cq = ea^d^ (with dual e“ = e“^dx^) and the 
connection uj\ = ui°’f^^dx^ = We mention that the torsion as dehned in (2.6) is 

a tensor under local Lorentz transformations and under diffeomorphisms. The vielbein is 
assumed to be orthonormal, which implies that the metric is given by 


9^iU — llab 6 6 t/ ) (^•'l) 

where rjab = diag(—1,1, ...1) is the Minkowski metric. The metric is assumed to be compat¬ 
ible with the connection namely r]ab\c = 0; thus Uabc = —^bac where | denotes co¬ 
variant differentiation with respect to The teleparallelism condition is imposed by as¬ 
suming that the curvature of 0 ;“^^ is zero, i.e. = 

0 (note however that the curvature constructed from the Christoffel connection is 


- 4 - 













not zero). One way to realize this condition is by considering the Weitzenbock connection, 
which is defined in terms of the vielbein in all coordinate frames as and 

therefore = 0 in the preferred frame Ca that defines the connection. 

Then, the following torsion scalar can be defined as 


T = (2.8) 

= , (2.9) 


where the tensor is defined by 

S^.ux ^ ^ \{g^^TpP'^ - gP’^TpP^) = ( 2 . 10 ) 


and 


TjyX^ T^y\) - ^fluX ^fluX 


( 2 . 11 ) 


is the contorsion tensor. In the teleparallel formulation, torsion is used to express the 
action (2.1) in an equivalent form (up to boundary terms) as 


Smin = -^^ j d^xeT- j d^xe(^^gP''df,4idy(j) + V^ . (2.12) 

Since T differs from R by boundary terms, namely T = —R + with ; denoting 

covariant differentiation with respect to the Christoffel connection, the above action (2.12) 
indeed gives exactly the same equations with the action (2.1). For example, the hairy black 
hole solution (2.2), (2.4), (2.5) can be obtained from the action (2.12) in four dimensions. 

In this work we are interested in finding local spherically symmetric solutions in which 
the presence of the torsion could leave an explicit signature on the solution, other than 
from the trivial one found from the action (2.12). To achieve this, the action (2.12) has to 
be modified in such a way as to include a coupling of torsion to the scalar field. 

The easiest way to modify the action (2.12) is to add a direct coupling of the scalar 
field to torsion 

Snon-min = J d^xcT- J d^x 6 df,(j)d^cj) + V + Cf {(t))T^ ■ (2.13) 

The resulting theory is different from its corresponding curvature counterpart and it leads 
to novel features when applied to a cosmological geometry [55-59]. However, it is not easy 
to apply the above theory to a spherically symmetric geometry in four dimensions. The 
reason is that for a diagonal vierbein the coupling of the scalar field to torsion gives a new 
(r — 0)-equation which acts as a strong constrain to the coupled field equations. Thus, the 
arising geometry is trivial and does not give hairy spherically symmetric solutions. Note, 
that in three dimensions the absence of this off-diagonal field equation allows hairy black 
hole solutions to be found [60]. 

Another possible extension of the action (2.12) is to include a derivative coupling of 
the scalar field to torsion. In particular, we consider such a coupling with the torsion scalar 
T, namely 


S = - 




J d^xeT — j d^xe 


,4)du4) + V 


(2.14) 
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The coupling parameter ^ has dimensions of length square, and therefore, -\/i^ introduces 
a new length scale in the theory. 

In order to extract the field equations, we need to perform variation with respect 
to the vielbein and the scalar field. Additionally, and more formally, in order to 
implement the teleparallelism condition Rabcd = Oj one adds in the action (2.14) the term 

Rated, where is a Lagrange multiplier. Thus, variation with respect to 

naturally provides Rated = 0, which can be considered as the equation of motion for the 
connection Variation of the action (2.14) with respect to the connection gives 

the equation of motion for the Lagrange multiplier. Since 



where 4”^ = Ignoring the boundary term and setting 5^3 = 0, we obtain 


+ vs;:) - (J—= 0 . (2.17) 


We mention that the above expressions contain also the connection ooate- Assuming for 
simplicity the Weitzenbdek connection, that is setting ujate = 0, which is a solution of the 
equation of motion Rated{‘^°‘te) — 0) acquire the field equations containing only the 
vielbein, namely 


2 

+ 4 ^ 




Additionally, variation of the action (2.14) with respect to the scalar field 4 gives 

6^3 = y'd^x|[e(l-2er)</.’^] _^-e^|<5</.- y'd^x[e(l - 2er)(/.’'^,5</-] , 

from where the scalar-field equation of motion is obtained as 


dV 


[<l-2Wl,p-e^ = 0. 


(2.18) 


(2.19) 


( 2 . 20 ) 
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Summarizing, equations (2.18) are the gravity field equations, while equation (2.20) is 
the Klein-Gordon one for the scalar field. Observe that T being a scalar can be identified 
with the Hubble parameter in a cosmological Friedmann-Robertson-Walker (FRW) back¬ 
ground as r = Goo = This implies that the derivative coupling on a 

FRW background will give the same cosmological evolution as the derivative coupling of 
the scalar field to Einstein tensor G^’^dfj^cpducj)- However, general cosmological perturba¬ 
tions will differentiate the two models because they will excite also the spatial components 
of the Einstein tensor. 

In the next section we will study local spherically symmetric solutions of the held 
equations (2.18) and (2.20) and we expect to hnd solutions with features different than the 
local black hole solutions of theories with a scalar held coupled to Einstein tensor with a 
derivative coupling [50-54]. 


3 Spherically symmetric solutions 


In this section we look for spherically symmetric solutions of the theory with a derivative 
coupling of the torsion to the scalar held in four dimensions. We consider the general 
spherically symmetric metric ansatz of the form 

ds'^ = -N{rfdt^ + K{r)-^dr‘^ + R{rfdn^ , (3.1) 

where + sin^O is the two-dimensional sphere, and N{r), K{r) and R{r) are 

three unknown functions. Note that we do not impose the relation N{r) = K{r) keeping 
the metric ansatz as general as possible. We could have chosen the gauge with R{r) = r, 
but we hnd it more convenient not to impose a particular gauge and to let also the function 
R{r) free. We consider the following diagonal vierbein which gives rise to the above metric 

= diag {N[r), K{r)~^, R{r), R{r) sin6») , (3.2) 


which is adequate to capture the main features of our theory. 

Inserting the vierbein (3.2) into the dehnition of the torsion scalar (2.8) we obtain 


T(r) = -2K‘^ 


R'rR' 2N' 

+ 


R \R 


N 


(3.3) 


where primes denote derivatives with respect to r and we denote k? = = SttGn. Also 

the held equations (2.18) become 


1 



1 




+ 2(2C(j)'^ - 



+ 2 ^^ 
R K 






R'/R' 2NG 1 


+ 8 & 


R'/R' 2N'\ 

'r\^ 


= 0 , 


(3.4) 

(3.5) 
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Equations (3.4)-(3.7) are accordingly the (A,/i) = tt,rr,99,9r components of the system 
(2.18), while the other components vanish. The Klein-Gordon equation of the scalar field 
(2.20) becomes 


I^KNR^ 


1 + A^K' 


.R'/R' 


R \R 



NR^ dV 
K 


(3.8) 


The off-diagonal equation (3.7) is a new feature of four or higher-dimensional telepar¬ 
allel gravity and it does not appear in the corresponding curvature-based equations of 
motion. It rises from the specific metric ansatz (3.1) and the vielbein (3.2) in spite that 
they are diagonal. Its appearance constrains the system of field equations making them 
more difficult to be solved. Its physical signihcance is that it connects directly the scalar 
field to the metric function K{r) of (3.1). 

Equations (3.4)-(3.8) are invariant under r-reparametrizations, in the sense that when 
r —)• f(r) and K —)• N^N,R^R, the equations remain form invariant. 

This implies that one of these equations is expected to form the constraint of the theory, 
as will indeed be shown in the following. More specifically we have four functions, namely 
N{r),K{r),(p(r),V{r) (the function R{r) is not considered in the enumeration because it 
corresponds to the choice of the radial gauge) and there are five equations of motion (3.4)- 
(3.8), but one of them is the constraint. So finally, we end up with four true equations 
for four unknowns. This is an interesting feature because the solution of the system for 
the diagonal vielbein will be basically unique without the need of any extra assumptions. 
On the contrary, in the curvature-based theories, for the same ansatz (3.1) the absence of 
the off-diagonal equation (3.7) provides three equations for the four unknowns. Therefore, 
one function has to be fixed arbitrarily, making the solutions quite sensitive to initial 
assumptions. 

Assuming a non-trivial profile for the scalar field, equation (3.7) is integrated to give 

= ( 3 . 9 ) 

where v is an integration constant with dimensions of inverse length square. This is ex¬ 
pected to characterize the charge of the scalar field. In section 4 we find an exact solution 
where the constant v is related to the derivative coupling ^ which determines the strength 
of the interaction of torsion to matter. However, in section 5, in the linearized solution 
found, these two parameters remain unrelated. 

Using (3.9) equations (3.4)-(3.6) and (3.8), along with the replacement of 
from (3.4) into (3.6), acquire the following form 

2R" A'2 1 

^ ^ + 2(1-2^k2j,2) ^2 (3.10) 
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R'/R' 2N'\ / 

^2 'rVr '^l^J ~Y 
N” R" N'/R' K'\ R'K' 

If ^ ^ ^ 7 j\'r^1{ ) ^ "r ^ 


NR^ 


1 + 4(K^ — 


R' /R 


R\R 


1-2^«;2z/2 

^2 


2(1-2^k2i/2) 



1 

R^ 


K2 


= 0 


2V+iy^ 

K2 


= 0, 


NR‘^ dV 
vK d(j) 


(3.11) 

(3.12) 

(3.13) 


The above system of equations is quite complicated and we aim, after performing a 
series of transformations to obtain a single, decoupled, master equation from where we can 
get the behaviour of the system. To do so, we consider as independent argument the variable 
(j) instead of the coordinate r and we find for any /(r) that f = -^Z, /" = ■^(/ —/^), 
where / = Then the system of field equations (3.10)-(3.13) becomes 


R R2 I k^{2V+u'^) 

R ^ 2R^ ~ 2f^ ^ 4z/2(1-2^k2z/2) “ 
R/R 2N\ / 1/2 1-2^k2j^2 1 \ ^ 

'r\R^ ^^ 1/2(1-6^k2j,2) y ~Y ^ R^J ~ 

N N R 1 k‘^{2V+i'‘^) 

N NR ~ 2R^ ^ 2v‘^R^ ^ 4z/2(l-2^At2z/2) “ 


f 2 

9 R 

{R 2n\ 

\ 

NR^ 


[-Nir)\ 

\ 


(3.14) 

(3.15) 

(3.16) 

(3.17) 


Replacing (3.16) by the sum of (3.14) and (3.16) and defining 


X = In R , 
R 


y = 


z = 


R ’ 

{RN^y 
RN^ ’ 


equations (3.14)-(3.17) become 


. 3 2 e-2^ /t2(2R+i/2) 

^~^2^ 2i/2 4 z/2(1-2^k2j,2^ 


= 0 , 


1 —2^k2z/2 _2^ 1/2 \ _ 

I/2(l-6eK2l/2)l ^2 e 

.. I 1,2 I 3k^{2V+u^) 

^ 2 ^ 2i/2 ^ 4 z/2(1-2Ck2i/2) 




1 


4^1/2 


Tya: =0. 


Thus, from (3.21), (3.23) we have 


(yz)' + 




k2(2R+i/2) 
2z/2(1-2^k2i/2)_ 


(3y+z) + ^""{y-z) = 0, 


(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 

(3.25) 
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while (3.22) gives 


K/ Kj 


(3.26) 


2 

Multiplying equation (3.26) by ~and adding this to equation (3.25) we get a 
new equation, in which when V is replaced from (3.22), we obtain equation (3.24). This 
reveals that equation (3.22) is a constraint, and thus, equation (3.24) is redundant. 

Note that the system of equations (3.4)-(3.6) for vanishing scalar field (j) = 0 and a 
cosmological constant as the potential V has the standard Schwarzschild-(A)dS solution. 
However, in the equivalent system (3.21)-(3.24) this limit has been lost, since setting v = 0 
in this system we get inconsistency. 

Defining the new variables 


Y = y^ (3.27) 

Z = z^, (3.28) 


and using the field equations (3.21)-(3.23), we get a differential equation for Y{x) as 


fY 

dx‘^ 


0. (3.29) 


Equation (3.29) consists the master equation which we were looking for, while the technical 
details for deriving this equation are given in the Appendix A. The quantities rj, fj are 
defined by equations (A.8), (A. 15). The strategy which we will follow is to find a solution 
of this equation and then use the other decoupled equations to specify the metric, the 
scalar field profile and the potential. 

In particular, finding the function Y from (3.29) allows us to use equations (A.4), 
(A. 16) in order to find Z as 


Z = 


4r?2y 



(3.30) 


Moreover, combining equations (A.10), (A.14), (A.16), we find the potential V as 


V = 


-2x 1 /dY 

-e - y - — ( — + 3y 


2 2r]\ dx 

Additionally, from equation (A.l) and the definitions (3.18), (3.19) we have 


(3.31) 



(3.32) 


which can be integrated to provide Then, we can use equations (3.20), (A.2) and 

(3.32) to obtain 


'dln{RN^)Y 

dx 


Z{x) 

y(x)’ 


(3.33) 
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the integration of which provides N{R). Finally, using (3.9), (3.18) and (3.32) we can find 

’ (3-34) 

and therefore the metric (3.1) acquires the form 

ds^ = -N^dt^ + -4^ + ■ (3-35) 

Although we have managed to extract from the initial coupled system of field equations 
the master equation (3.29) plus decoupled equations, the general solution of (3.29) is still 
difficult to be determined analytically. In the next two sections we will give two solutions of 
the theory and discuss their physical significance. The first one is a special exact solution 
of the field equations where the potential turns out to have a simple form. The second one 
is obtained from an asymptotic solution of the master equation (3.29). 


4 A new class of wormhole-like solutions 

In this section we will find a special exact solution of the above system of equations, and 


in particular written in the form (A.3)-(A.5). Consider the potential 

V{x) = ae-^^ + (3 , (4.1) 

where a and /3 are constants. Then, from equations (A.6) and (A.7) we obtain 

Y(x) = cie"^^ + - |(i/^ + 2/3) , (4.2) 

Z{x) = C 2 e“^ - - 6i]c4j 6“^"^ - 377 ( 1 /^ + 2/3), (4.3) 

where ci, C 2 are integration constants. Choosing ci = C 2 = 0 we find 

Y{x) = Ae-2^ - B , (4.4) 

Z{x) = Ce-^^-D, (4.5) 


where 

^ 1 

A = ^ - 27ja , 
il = |(i/2 + 2/3) , 

C = Gua- K , 

D = 9B. (4.6) 


Thus, inserting (4.4), (4.5) in (A.4), we find that this equation is satisfied by choosing 


Case I : a = 
A = 


5 
8^2 
4k2^ 




2k^ 


1 


8k2^ ’ 


B = — , C = , D = Gk^ , 

3 


(4.7) 
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or 


Case II 


- 1 «_ 2_ 1 

“ 2^2 ’ ^ "" 5^2^ ’ 


A = 


6 


, B = 


5^2 

~6 


, c/ = 




D = 


15«;2 


(4.8) 


Since y, Z > 0, it is implied that Ae~‘^^ — B > {), Ce“2* — D > Q, and therefore 
d., C > 0, which requires ^ > 0. Observe that in both cases, because the constant (5 of the 
potential (4.1) is proportional to 1/^, it is scaled as and therefore acts as a positive 
cosmological constant. Also note that the charge of the scalar field v is proportional to 
Moreover, we deduce that R < y^, \J~%- This means that for the case (4.7) it is 
R < \/2^5 while for the case (4.8) it is ii < v^. 

Now, knowing the solution for Y{x) from (4.4), we can use equation (3.32) to find the 
scalar field as 


ei 


(j){R) = 4>i + —j= arctan 


A 1 


- 1 


(4.9) 


where ei = ±1 and (pi is an integration constant. It proves convenient to redefine the 
held 0, absorbing the integration constant cpi and the index ei. In particular, setting 
(p{R) = ei[(p{R) — (pi] we result to 


1 


(p{R) = —!= arctan 




A 1 


- 1 


(4.10) 


Thus, for the case I it is 

kR) = 

while for the case II we have 

kR) = 


3 1 

2 K 


arctan 


i?2 


- 1 


(4.11) 


6 1 
5 K 


arctan 


i?2 


- 1 


(4.12) 


The potential V can be now reconstructed by inverting (4.9), namely 

\Vb{(P-(Pi) 


^(<(’) = ^ tan" 


aB 

+ y^ + /3, 


or 


= ^tan^ (v^(^) + ^ + /^- 


Hence, for the case I we have 




K(p + 


16^«;2 ’ 


(4.13) 

(4.14) 

(4.15) 
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while for the case II it is 


V{^) = 




■ tan 



(4.16) 


The dependence of F to i? for the case I is 


1 




8^ k 2 8^2 ^2 ’ 


(4.17) 


while for the case II 


1 


yiR) = 7772 + 


1 


(4.18) 


5^k2 2k? W' 

Note that the solutions for the scalar field and the potential have analogous forms with 
the corresponding solutions of the scalar field and potential of the MTZ solutions (relations 
(2.4) and (2.5) respectively). However, while the cosmological constant A = K^V{<p = 0) of 
the MTZ solution is an arbitrary parameter, here, the corresponding cosmological constant 
is given for case I by A = and for case II by A = Thus, the gravity-scalar held 
equations have hxed this constant to depend on the coupling of torsion to the scalar held. 
Not only this, but also the integration constant v of the scalar held is determined from the 
coupling ^ through Eqs. (4.7), (4.8). This feature also occurs to the solutions of scalar- 
tensor theories, where usually the scalar charge of a local hairy black hole solution is related 
to the other scales of the theory like mass or charge, unless a symmetry is spontaneously 
broken, like conformal symmetry [38]. 

The lapse function can be calculated from equation (5.11). Using the parameter values 
from (4.6), (4.7), (4.8), we hnd 

d\n{RN‘^) 


dx 


= 3e. 


which has as a solution 


N[R) ^ 

with e = ±1 and C an integration constant. Therefore, 


N = CR for e = 1 


or 


N = 




for e = — 1. 


Finally, inserting the above relations into (3.35), we hnd the metric 

dR^ 


ds^ = -C'^R^^-^dt^ + 


u'^{A-BR^) 


+ R^dn^ 


e = ±1. 


(4.19) 


(4.20) 


(4.21) 


(4.22) 


Thus, for the case I it is 


ds^ = -C^R^^-^dt-^ + + RW 

-tX/ 


(4.23) 


while for the case II we have 


ds‘^ = 


, 2 6.^ dR^ ^2 7r-.2 

dt + + R dn^ 




(4.24) 
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These metrics, together with (l){R) from (4.10) and V{(p) from (4.14) form the exact 
spherically symmetric solutions that we have found. The solutions contain the two cou¬ 
plings K, ^ and one integration constant C. Note that when the metric function in (4.23) or 
(4.24) has a singularity, i.e. at the maximum value = 2^ or R^ = the accompanied 
scalar field (p vanishes, = 0. Moreover, at the origin R = 0, the scalar field is still hnite. 
Therefore, the scalar field is regular everywhere. However, the potential V becomes infinite 
at the origin i? = 0, as it is seen from (4.17), (4.18). Finally, the Ricci and the Kretchmann 
scalars diverge at the origin but are finite elsewhere. For example, for the metric (4.24) 
and e = 1 it is R^ukxR^’'^^ = ■ 

Defining b{R) = v^BR^ -|- (1 — v‘^A)R the metric (4.22) can be written in the form 

ds^ = -C^R^^-^dt^ + ^ + R^dQ^ , (4.25) 

where for Case I it is h[R) = ^R^ + §7?, while for Case II it is b{R) = ^7?^ + |7?. The 
lapse function N{R) is always regular, therefore, the metric (4.25) does not represent a 
black hole. It has the form of a wormhole [61, 62] with “throat” at 77 = AjB. The 
difference is, however, that here the solution extends inside the “throat” up to the origin 
77 = 0, while in a usual wormhole the solution extends outside the “throat” up to the 
“mouth” or up to infinity. The behaviour of this wormhole-like solution is controlled by 
the length scale ^ and this solution merits more investigation. 


5 Spherically symmetric solutions at large distances 


In this section we will solve the master equation (3.29) perturbatively. We will assume that 
for X —)• oo, where x is related to the radial distance as (3.18), the function Y{x) has the 
expansion Y{x) = Yq + Yi{x) -|- ..., where |Fi(x)| <C |lo|- So, we will obtain large distance 
linearized solutions. 

At zero order (3.29) gives 


>0 = 


3(2o--l) ’ 


(5.1) 


where we have introduced the constant 



V 


(5.2) 


with e = ±1. Since, due to (A.l), it is T > 0, we have To > 0 and therefore r?(2cr — 1) > 0. 
Obviously the exponential factor in (3.29) is suppressed compared to the asymptotic value 
Yq. To continue with the linearized correction Yi we assume additionally that |li(x)| <C 
[To + ^^1- The master equation (3.29) at first order, after substituting the zeroth order 
value (5.1), gives 


+ 


4(l-a) + (l-2u) 


o— 2x 


^ + 3(1-2a) 
dx 


l + (l-2a) 


o— 2x 


Yi + 


2 ( 76 - 2 ^ 

3i/2 


= 0 . 

(5.3) 
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Assuming further that R S> ^ i-a\ ’ equation (5.3) becomes 


(fYi ^dYi , 2cje-2^ 

^ + 3(l-2a)Fi + ^ 


= 0 . 


The general solution of (5.4) is given by 

Ti = CAe^" + + 

where cx , are integration constants and 


2a 


„-2x 


3zy2(l-2f7) 


A = 2(u-1)-a/1-2u+4o- 2 < 0 , ;ti = 2(o--l) +V l-2a+4o-2 


(5.4) 


(5.5) 


(5.6) 


with // > A being either positive or negative. Depending on the possible parameter values 
we have two cases: 

• Case A 

For fj < ^ we have rj < 0. Then the only consistent case is for e = —1, ^ 

(thus the coupling ^ must be positive), and finally a < — |. The values of A,/r are 
A < —2, —3/2 < ^ < 0, and therefore the total linearized solution becomes 

, 2a 


Y = 


3(l-2a) 


+ + 


3i^2(i_2o-) 


e-2^ + 


(5.7) 


with 


0 < m = 2(1 —a) — x/l — 2a+Aa‘^ < - 
I = 2(l-o-) + \/l-2o-+4cr2 > 2, 


(5.8) 


or in terms of the radial distance R 
2\r]i\v‘^ 


Y = 


+ CmR + 


2a 


-R-^ + ciR-^, 


(5.9) 


3(1-2a) ' ' 3zy2(l-2a)‘ 

where C£,Cm are integration constants. According to the previous assumptions of 


linearization, the above solution is valid provided that 


Cm.R 


2—m 


2a 


_L Cl 

3u'^(l-2a) 


< 


3 T- 2 o-) '^^’ 3^-20 ^^' inequalities mean that R should be larger than some 

minimum distance, which decreases as Cm decreases. 

Using the above solution of the master equation and following the procedure described 
in Section 3, the metric (3.35) is found to be 

dR^ 


ds"^ = -N^dt^ + 


2a 


r?2 I ^ f>2-m I __ 

3(1-20-)-^ 3(1-20-) 


I 

' dP— 


+ R^dn'^ , 


(5.10) 


where are redehned integration constants. The lapse metric 

function N can be obtained from the equations (3.30) and (3.33) from which we 
obtain 

t iRV ..N 

(5.11) 


dln(AA^2) ^ 
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where C = ±1 is a sign symbol. Integrating equation (5.11) we extract the lapse 
function as 


where 


J{R) = I 


dR 


N^{R)= 

R L 


du 


, A. 


RY{R) 




1_7n 

-U + CmU 2 + 


2a 


3(1-2(t) 


+ 


ce 


(5.12) 


(5.13) 




3(1-2(t) 

and c > 0 is an integration constant. 

The potential can be determined after the insertion of the solution (5.9) into (3.31) 


to find 


V = 


(l+2a)iy^ 


3—m Cr 


3-8a 


o-ou 1 3-£ Cl 

2(1 —2(t) 2r\v‘^ RY- 6(1 —2(7)r/z/2 2'qiy‘^ R^ 

Finally the scalar field conhguration can be determined inserting the solution (5.9) 
into (3.32) to obtain 

dR 


= </)! + eiZ 2 


/ 


3(i-2o-) 


i? 2 + CmR^-'^ + 


2a 


__£L_ 

3(1-2(t) 


(5.15) 


where ci is another ±1 sign and 4>i is an integration constant, asymptotic Going 
back to the potential (5.14), we can find the next-to-leading order correction of V as 
a function of the scalar field, namely 


V(6) ~ + _ (3 m)Cm 

^ 2(1-2a) 2r/z/2 


(5.16) 


where (^ = ei ((/> — cj)i^. 

The integral (5.13) cannot be performed analytically, however, the leading asymp¬ 
totic behaviour is found to be J{R) ~ Ini?, therefore, the leading asymptotic 

behaviour of the lapse function is 


N‘^{R) 


2(7 


2|7/|z/2 

3(l-2o-)J i?i-3C • 


(5.17) 


Of course, the exact integral (5.13) gives more information at even smaller distances. 
The asymptotic form of the metric (3.35) becomes 


dsoo = -c 


1 


2|?7|i22 

3(l-2a)J i?i-3< 


dR + 


3(1-2c 7) di?2 
2|r?|i/4 W 


+ R^dn^. 


(5.18) 


Therefore, for the branch with C = 1) rescaling the time t to t = ^ [ 3 ^ 1 - 20 -) ] ^ 

we find the asymptotic metric 


dst = -,^^!^R^dP + + R^dn^ . 


3(l-2a)' 


2\r]\u^ R? 


(5.19) 
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The resulted metric (5.19) indicates that the space is asymptotically AdS and some 
interesting features appear. The length scale introduced here is given by = 
3^1-20-) corresponding effective cosmological constant is Aeg = —3/^. On 

the other hand, the potential V of the solution asymptotically has a different negative 
constant given by Eq. (5.14) for R —>• oo. Therefore, the constant coming from the 
potential of the action and the effective cosmological constant of the asymptotically 
AdS spacetime can have different values, and this is due to the fact that the scalar field 
acquires a non-trivial profile at infinity. In this way it introduces another energy scale 
asymptotically. Note that the ratio of Agfr and the constant appearing in the potential 
is proportional to therefore, if ~ it happens in the exact solution 

we found, the two constants are of the same order. This behaviour is different from 
the exact MTZ solution where the potential asymptotically becomes a cosmological 
constant which is the same with the cosmological constant of the asymptotically AdS 
spacetime seen in Eqs. (2.2), (2.3). Hence, in our case, the presence of the scalar 
field and its coupling to torsion modifies the asymptotic form of the spacetime. 

In the linearized large distance solutions we discuss in this Section the length scales, 
defined by the cosmological constants discussed above, depend on two quantities, 
the derivative coupling constant ^ and the integration constant u which expresses 
the scalar charge. Thus, contrary to the previous wormhole solution, in this case 
the two constants are not interrelated and this is explained since the linearized so¬ 
lution found is a general asymptotic solution with the correct number of integration 
constants. Note also that the derivative coupling constant ^ enters the solution only 
through the specific combination while the integration constant u appears also 

independently. This combination is a remnant of the interaction term 
in the initial Lagrangian. Therefore, the value of ^ is not significant on its own, but 
only in combination with the integration constant v of the scalar field. 

We note here that the aforementioned AdS behaviour, as opposed to the dS behaviour, 
creates an (extra) attraction force at local level at large distances, which may be useful 
at galactic or cluster of galaxies scales, and in particular can offer a novel mechanism 
for their stability. Also in the metric (5.10) the term 1?^“™ expresses another large 
distance scale where gravity is modified, and for instance, in the case where m = 1 
this term becomes proportional to R. Then in [63] it was claimed that in Weyl gravity 
a linear term R can successfully explain the galactic rotation curves without the need 
of adding dark matter. 

• Case B 

Eor u > I we have g > Q. Then the only consistent case is for e = —1, 

2 

(thus the coupling ^ must be positive) and The values of are // > 0, 

—2 < A < —3/2, and thus it should be = 0 otherwise Yi keeps increasing at large 
distances and the condition |li| ^ [ToI is not satisfied. Therefore, for u > ^ the 
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solution becomes 


Y = 


2rju‘^ 


3(2c7-1) 


+ - 


2(7 


„-2x 


3i^2(2ct-1) 


with 


or in terms of R 


- < ^ = 2(l-o-) + Vl-2a+4o-2 <2, 


Y = 


2rju^ 


+ ciR-’^ - 


2a 


-R 


-2 


3(2cr-l) ' 2,u'^{2a-iy 

where q is an integration constant. The metric (3.35) takes the form 


(5.20) 

(5.21) 

(5.22) 


ds^ = -N^dt^ + 


dR^ 


2r]u‘^ 


3(2(t-1) 


i?2 + CiR^-^ - 


2(7 


+ R^dn^ 


(5.23) 


3(2(7-1) 


where C£ = z/2q is a redefined integration constant. Following the same procedure 
as in in case A we find similar forms for the lapse function and the scalar field as in 
case A, while the form of the potential is found to be 


(1 + 2cj)i^2 3 — 771 Cm 3 —8cJ 1 3 — i Cl 

2{l — 2a) 27/1/2 j^m 6(1 —2f7)r/z/2 i ?2 2r/z/2 R^ 


The physical significance of this solution is similar to the discussed case A solution. 


6 Conclusions 

In this work we investigated spherically symmetric solutions in scalar-torsion gravity theo¬ 
ries. Using the teleparallel formulation of GR, which presents the torsion scalar field T, we 
allowed a derivative coupling of a scalar field 4> to the torsion scalar. The resulted theory 
is a novel scalar-torsion theory different from the corresponding curvature-based gravity 
theories. It has the advantage of giving second order field equations leading to a healthy 
theory without ghosts. 

After extracting the general field equations, we applied them in the background of a 
spherically symmetric geometry. A novel feature of the arising coupled system of gravity- 
scalar field equations is the appearance of an off-diagonal equation, which couples the scalar 
field directly to a metric function. The appearance of this equation constrains further the 
system of the coupled held equations, making them difficult to be solved analytically. In 
spite of that, we managed to extract a single decoupled master equation, the solution of 
which allows to hnd one by one all the other unknown functions. 

Although the solution of the system of equations is unique for the metric components, 
the scalar held prohle, and the scalar-held potential, we have found a specihc form for the 
scalar potential which provides an exact solution of the held equations. This solution is 
special from the point of view of the number of integration constants, and as a result, both 
the integration constant of the scalar held and the asymptotic value of the potential depend 
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on the derivative coupling constant. Interestingly, the obtained solution is a wormhole-like 
one, however it extends inside the “throat” of the wormhole up to the origin, while in 
a usual wormhole the solution extends outside the “throat” up to the “mouth” or up to 
infinity. Moreover, the accompanied scalar field remains regular both at the origin and at 
the “throat”. Hence, this solution belongs to a new class of wormhole-like solutions which 
can have interesting physical features. 

Using the master equation, we obtained linearized solutions that are valid at large 
radial distances. Interestingly enough, we showed that at large distances the space is AdS, 
which implies that in the model at hand, at large distances we obtain a negative effective 
cosmological constant. However, the asymptotic form of the potential provides a different 
negative constant. This is a new feature, not present for example in the MTZ solution, 
and is due to the derivative coupling along with the non-vanishing asymptotic scalar held. 
Hence, in our case, the presence of the scalar held and its coupling to torsion modihes the 
asymptotic form of the spacetime. In these solutions the derivative coupling constant and 
the integration constant of the scalar held are not interrelated, contrary to the previous 
wormhole solution, and this is due to the fact that the solutions are general from the point 
of view of the number of integration constants. 

The fact that the derivative coupling creates an anti de-Sitter space in the present 
solution, instead of a de-Sitter one, gives an attractive term around a spherical object, 
which might have interesting implications on structure formation and cosmology. To the 
same direction may act the fact that apart from the above behavior, the radial metric 
potential acquires another contribution at smaller distances that can be close to linear-in- 
distance, and thus playing the role of an effective dark matter component. 

In summary, as we observe, the derivative coupling in the framework of teleparallel 
gravity brings novel features that cannot be obtained in the usual curvature nonminimal 
couplings. However, a crucial test for the scalar-torsion theories to give distinct features 
from the scalar-tensor theories of GR, is to find exact solutions representing a central 
object, like a black hole with scalar hair. Before considering this theory as a candidate for 
the description of nature it has to be shown that it is viable phenomenologically and the 
model parameters have to be determined. 

A Calculation of the master equation 

In this appendix we will give the technical details for obtaining the master equation (3.29). 
Defining the variables 


Y = y^ (A.l) 

Z = , (A.2) 
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we find that ij = and therefore equations (3.21)-(3.23) become 


dY e 

—+ 3y- 
dx 

dZ 


-2* n^{2V+-u'^) 

+ ^ n n. = 0 




YZ = 


^-2x 


dx 


+ Z- 


3K^{2V+iy^) 
z /2 ' 2 i ''^{1 — 2 ^ k '^ u '^) 


+ 


= 0 . 


Equations (A.3), (A.5) are linear differential equations, whose solutions are 
Y{x) = e-3^ 


Cl + ^ - ^z^2g3x - 2r/ y Ve^^dx 


^—x 


Z(x} = e ^ ^C 2 - 2 -3r/i/2e^ ~ J ’ 

where ci, C 2 are integration constants and we have defined the constant 

V 


2i^2(l — 2^F^u'^) 


(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 


Finally, inserting (A.6), (A.7) in (A.4) we obtain an integro-differential equation for the 
potential V, considered as a function of x, i.e. V = V{x). Therefore, if the potential 
V{x) was given, in general the system would be inconsistent, and the reason behind this is 
the extra off-diagonal equation (3.7). However, if the potential V{x) satisfies this integro- 
differential equation the system is consistent. 

Let us transform the aforementioned integro-differential equation into a differential 
equation. Dehning the function 




we immediately obtain 


E = e 


— 


dh 

dx 


(A.9) 


(A.IO) 


while fVe^dx = fe "^^j^dx = e ‘^^h + 2 f e “^^hdx. Using /i(x), equations (A.6) and 
(A.7) give 

Y{x) = (ci + J - 2r?h) , (A.ll) 

Z{x) = e~^ ^C 2 — ^“2 3r/z/2e^ — 6rie~'^^h — I2r] Je~‘^^hdx'^ . (A.12) 


Thus, replacing Z from (A.4) in (A. 12) and differentiating we get 


d 


dx 


.F2 


Y 


_(e^—) = _ -3r?i/2e^ 




dh 


dx 


Additionally, from (3.22) we deduce that 


F = h[ 


( 1 
2u‘^r]^ 


2x 




—Sx 


dh 

dx 


(A.13) 


(A.14) 
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where we have introduced the constant 


Therefore, from (A.11), (A.14), the function F can be expressed as 




while (A. 13) becomes 




dx\ Y 


dx 


+ 9e^Y - —e- 


(A.15) 


(A.16) 


(A.17) 


Hence, given (A.16), equation (A.17) is a second order autonomous equation for Y (x) which 
is also written as 


2F 




dx \Y 


dx 


+ F^ - 9Y^ + -^e-^^Y = 0 . 


(A.18) 


Finally, using (A.16) and (A.18) we extract the differential equation for Y(x} as 




dx'^ Y \ dx J 


dY 




2 ^-2x 


, , ■+ 3Y+ 27112^-122-Y- _ 

Y J dx rj^ ^ _g 3 y + 2 r]i 2 ^ 


= 0. (A.19) 


In summary, we have managed to replace the complicated system of field equations 
(3.4)-(3.8) with just the equation (A.19) plus decoupled equations. 
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